Role of interactions in the far-infrared spectrum of a lateral quantum dot molecule 



o 
o 

(N 

Oh' 
(D 

m 
(N 



X3 

o 



> 

a\ 
o 
m 
o 



I 

o 
o 



X 
S3 



M. Mario, A. Harju, and R. M. Nieminen 
Laboratory of Physics, Helsinki University of Technology, P. O. Box 1100 FIN-02015 HUT, Finland 

(Dated: February 2, 2008) 

We study the effects of electron-electron correlations and confinement potential on the far-infrared 
spectrum of a lateral two-electron quantum dot molecule by exact diagonalization. The calculated 
spectra directly reflect the lowered symmetry of the external confinement potential. Surprisingly, we 
find interactions to drive the spectrum towards that of a high-symmetry parabolic quantum dot. We 
conclude that far-infrared spectroscopy is suitable for probing effective confinement of the electrons 
in a quantum dot system, even if interaction effects cannot be resolved in a direct fashion. 
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Nanoscale semiconductor structures are very promis- 
ing for future components of microelectronic devices. 
These systems are also scientifically very interesting as 
they exhibit novel and fundamental quantum effects. The 
most prominent difference between the two-dimensional 
artificial atoms, or quantum dots (QD), and their normal 
counterparts are the enhanced correlation and magnetic 
field effects. In addition to that, the new features can 
be controlled by the tunable system parameters, both 
experimentally and in theoretical models. Coupling to- 
gether QDs, one can construct QD molecules (QDM). 
The electronic structure of these is also quite intriguing. 
For example, the two-electron QDM has a highly non- 
trivial spin-phase diagram and composite-particle struc- 
ture of the wave function 1]. One experimental possi- 
bility to probe the new physics of nanostructures is use 
far-infrared (FIR) absorption spectroscopy 0, 1^. For 
highly symmetric QDs, the FIR spectra reflect only the 
center-of-mass properties of the system. If the symme- 
try is lowered, the relative motion of the electrons starts 
to couple to the center-of-mass motion, enabling one to 
study the many-body correlation effects using FIR. 

In this Letter we present FIR spectra for a two-electron 
QDM using exact diagonalization. The method used 
solves accurately the quantum-mechanical model for a 
QDM, revealing the important electron-electron correla- 
tion effects. Our calculated spectra show clear deviations 
from the spectra of highly symmetric QDs. To analyze 
these deviations, we vary our system parameters in order 
to separate the correlation effects from the single-particle 
effects of the low-symmetry confinement potential. This 
analysis shows, however, that these effects are ultimately 
entangled. Furthermore, we find, very surprisingly, that 
the correlation effects actually compensate the lowered 
symmetry of the QD confinement potential, resulting in 
spectra much closer to an electron in a high-symmetry 
QD than in a low-symmetry one. 

The starting point for understanding FIR spectra is 
the highly symmetric, parabolic QD. This spectrum 
consists of two Kohn modes, whose dispersion does 
not depend on the number of confined electrons nor 
their interactions "^l. Several experiments on QDs have 



shown deviations of the FIR spectra from the Kohn 
modes 0, 0, S Si- It is clear that these deviations 
require a non-parabolic QD, but the detailed cause of 
the deviations, and thus the interpretation of the mea- 
sured spectra, is difhcult to obtain. An especially inter- 
esting question is how the electron-electron interaction 
and correlation effects emerge in the FIR spectrum when 
the symmetry of the QD is lowered. In our model, we 
can tune both the deviation of the potential from being 
perfectly parabolic as well as the electron-electron inter- 
action strength, which enables us to analyze the source 
of the non-trivial details seen in the FIR spectra. This 
type of analysis is extremely important for explaining the 
details of the spectra. 

The calculated FIR spectra have shown that a bro- 
ken xy-symmetry {e.g. elliptic or rectangular) is re- 
quired to observe the zero-field splitting of the Kohn 
modes as well as additional modes, and a broken ro- 
tational symmetry (e.g. sq uare) is required to observe 
anticrossings 0, 0, 0, 0, However, the role 
of interaction effects and lowered QD symmetry in FIR 
spectra is not fully analyzed in these studies. Our QDM 
model suits perfectly for this detailed analysis as the FIR 
spectra show all the above-mentioned features. One in- 
teresting finding is that the electron-electron interactions 
drive the FIR spectrum closer to the Kohn modes instead 
of making the deviations more pronounced. This means 
that the electrons feel an effective potential that is more 
parabolic than the the bare confinement one. Thus the 
interactions smoothen the non-parabolicity of the exter- 
nal potential. 

We model the two-electron QDM by a 2D Hamiltonian 
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where Vc is the external confining potential, for which 
we choose ^m*u}Q min{(a; - d/2)^ -I- y^, (x + d/2)^ + y^}. 
This potential separates to two QDs at large inter-dot 
distances d, and with d = it simplifies to one parabolic 
QD. We use the GaAs material parameters m* /me — 
0.067 and e = 12.4, and the confinement strength hwa = 
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3.0 meV. A is the vector potential of the magnetic field 
(along the z axis) taken in the symmetric gauge. We use 
C to scale the Coulomb interaction strength. The wave 
functions and energies of the ground and excited states 
are found by the exact diagonalization of the Hamiltonian 
matrix as in Ref. 0. 

Using the Fermi golden rule within the electric-dipole 
approximation for the perturbing electromagnetic field, 
the transition probability from the ground state to the 
rth excited state can be calculated as 
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where we have used right- or left-handed circular polar- 
ization of the field, corresponding to the two possible 
signs marked with '±'. We mainly present results for 
unpolarized light, averaging over the two circular polar- 
izations. In addition to that, we decompose one selected 
spectrum to two circular as well as to two axial polariza- 
tions. 

The calculated FIR spectra are shown in Fig.^for sev- 
eral different system parameters. We plot the transition 
energy as a function of the magnetic field using a line 
whose width is proportional to the transition probabil- 
ity. We have only included the transitions which have 
a probability of more than one percent of the maximum 
value. To make it easier to compare the relative strengths 
of the different resonances, we have plotted the probabil- 
ities also separately below each spectrum. In addition to 
that, we show the energies of the two Kohn modes of a 
parabolic QD using open circles. The upper row of figures 
present spectra for the total spin S = and the lower one 
for S = 1. The spectra for both spin types are shown for 
two inter-dot distances d and three Coulomb strengths 
C. The ground state of our system changes from S = 
to S* = 1 when the magnetic field is increased (see Ref. 
for more details), and this transition is shown by a ver- 
tical line in the four spectra of the fully interacting cases 
C = 1 . One should note that the transition probabilities 
between different spins are zero, and thus S is conserved 
in the allowed transitions. 

As a general feature of the calculated spectra shown 
in Fig. n one can see that each spectrum has as a major 
component two branches, where the higher one {uJ+) has 
a positive dispersion and the lower (lo-) a negative one. 
For the parabolic QD, these branches are the two Kohn 
modes (circles in the figures). Now the lowered symme- 
try of the QD and the correlation effects show up in the 
FIR spectra as a deviation from the Kohn modes. These 
deviations are the splitting of lu^ and lu^ at _B = 0, the 
anticrossings in and an additional mode a;+2 above 
a;+ in the 5=1 spectra. Next we will analyze these 
deviations in detail. 

Figs, n (a) and (b) show the (C = 1) spectra for the 
distance d = 10 nm which corresponds to closely coupled 



QDs. One can see that there are only minor deviations 
from the Kohn modes. For S = 0, one can see a small 
B = splitting of the two modes. There are also anti- 
crossings in the upper modes, but these are hardly visible 
as the energy gap is much smaller than the line width. 
These anticrossings are, however, visible as oscillations 
in the transition probabilities shown below the spectra. 
These probabilities also show that the upper mode of the 
S — I spectra actually consists of two levels. These levels 
are energetically nearly degenerate, but the probabilities 
vary as can be seen in lower part of Fig. ^ (b) . The sum 
of the two upper modes adds up to a probability which is 
nearly constant and close to the probability of the lower 
mode. 

If we increase the distance d between QDs, the confine- 
ment potential becomes less parabolic. In the extreme 
limit, however, the two QDs decouple and one is left with 
two electrons localized in separated dots. Figs.n(c) and 
(d) show the spectra for d — 30 nm which is in the in- 
teresting range of couphng between the two QDs. One 
can see that the features discussed for d = 10 nm are 
again present, but at this time much clearer in the spec- 
tra. The B = splitting is rather large for S* = 0, and 
there is also a small splitting in the 5* = 1 case. There 
are many more anticrossings than in the d = 10 nm case, 
and those for weak B now have a clear energy gap. The 
two nearly degenerate modes of the S — 1 case are now 
split in energy, and the higher one (w+2) has a smaller 
transition probability than 07+. The comparison of the 
spectra of Figs, n (a) to (c) and (b) to (d) shows that the 
first anticrossing point is rather insensitive to the change 
of d even if the gap in energy opens up clearly. A further 
comparison shows that at B = 0, the uj+ mode stays at 
the energy of 3 meV, but lu- is lowered. The splitting 
of the modes results from the fact that the symmetry in 
X- and y-directions is broken in the confinement poten- 
tial for a non-zero d. For our Hamiltonian, the confine- 
ment along y-direction is kept parabolic even for d > 0, 
and the strength in this direction is 3 meV. For the x- 
dircction, the potential is non-parabolic for d > 0, and 
the effective strength of an approximating parabolic con- 
finement is smaller than 3 meV for any finite d. In this 
way, the B = modes correspond to excitations to ei- 
ther the short or long axis of the system. This fact can 
be further demonstrated by considering the spectra of x- 
and y-polarized radiation, presented in Fig. [21(a). One 
can see that at _B = the transitions are given purely by 
either of the two linear polarizations. For a non-zero B, 
the two linear polarizations mix due to a term A • V in 
the Hamiltonian. Measurements using linearly polarized 
radiation lead to similar conclusions, see Refs. 0,[3,I3- 
For a strong magnetic field, the magnetic field length 
{Ib = ^ft/ei?) becomes smaller than the confinement 
length (Zo = \/fi'/^o'm*) at i? « 1.7 T. Thus at greater 
_B, the harmonic Landau potential dominates over the 
QDM potential and the electrons start to localize in the 
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FIG. 1: FIR spectra of QDM with two different distances [(a) and (b): d = 10 nm, (e)-(h): d — 30 ran] between QDs and 
full [(a)-(d)] and reduced [(e)-(h)] Coulomb strengths C. Upper subfigures correspond to S = spectra and lower 5 = 1 one. 
Vertical lines indicates the transition points from spin-singlet {S = 0) to spin-triplet {S — 1) state. Spectra show the energy of 
absorbed light as a function of magnetic field, and the width of lines indicate the transition probabilities, also plotted separately 
below each spectrum. Circles show the Kohn modes of an isolated parabolic QD. 
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FIG. 2: Polarization dependence of FIR absorption for a 
QDM with d = 30 nm, S = 0, and C = 1. The left panel 
corresponds to linearly polarized radiation and the right panel 
to two circular polarizations. At _B = the modes are linearly 
polarized and for large B they approach circular polarizations. 

minima of the confinement potential ij. Therefore, the 
non-parabolicity is diminished. The same effect can also 
be seen in the absorption of circularly polarized radia- 



tion, shown in Fig. |5] 

To separate the effects of the non-parabolic confine- 
ment potential from those of interactions on the FIR 
spectra, we have performed calculations with reduced 
Coulomb interaction strengths < C < 1. We present 
results for the d — 30 nm case as it shows highly non- 
trivial spectra. Figs.n(g) and (h) show a non- interacting 
spectra C — 0, and Fig. (e) shows a 5 = spectrum 
for C = 0.2 and Fig. [H (f) the S* = 1 one for C = 0.5. 
Different intermediate values of C are chosen for differ- 
ent S as the ones shown are the most representative ones. 
Comparison of the C < 1 spectra to those of Figs. ^ (c) 
and (d) with C = 1 shows a very surprising feature: the 
spectra for the reduced interaction strengths differ much 
more from the Kohn modes than the fully interacting 
cases. The splitting of the modes a.t B — grows as 
the interactions are reduced, the anticrossing gaps are in 
general clearer, although a direct comparison is not very 
easy as there are many shifts in the anticrossing posi- 
tions. In addition, the a;+2 mode of 5 = 1 becomes more 
distinct in the spectra as C is reduced. 

We start the analysis of the C < 1 spectra from the 
zero field splitting of the two modes. The potential along 
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the y-axis is parabolic, and tj+ at B = is not affected 
by the interactions. On the other hand, the excitation 
energy along the non-parabolic axis is clearly influenced 
by scaling C. The Coulomb repulsion effectively steep- 
ens the confinement which leads to an increase in the 
excitation energy. 

In experiments for elliptic QD lattices, the zero field 
resonance of uj- was found to exhibit an unexplained 
and dear dependence on the number of electrons in the 
QD 8] . The lowest energy was observed with the smallest 
number of electrons. As the non-interacting S = spec- 
trum coincides with the spectrum of a one-electron QDM, 
one can study the particle number dependence of uj- for 
the smallest particle numbers. Interestingly, we find a 
very similar dependence in our data by comparing the 
corresponding spectra in Figs, ^(g) and (c). The reduc- 
tion of the B = excitation in Ref. occurs only in the 
long-axis direction. This implies a non-parabolic poten- 
tial in that direction, and on the other hand, a parabolic 
one in the short-axis direction. 

Another interesting question is how the anticrossings 
change when the interactions are scaled down with C < 
1. For all interaction strengths, the energy gap at the 
anticrossing point gets smaller as one moves to stronger 
B. This can be seen to result from the localizing effect 
of the magnetic field: as the electrons localize around 
the potential minima, they feel a more parabolic exter- 
nal confinement. Now when the interactions are made 
weaker, the positions of the anticrossings move, in gen- 
eral, to higher B. However, the energy gaps associated 
to fully and non-interacting cases are nearly the same; 
the change is typically less than 10%. Thus the reduced 
interaction is compensated by the stronger localizing ef- 
fect of B at the new anticrossing position of higher B. 
This can be understood by noting that interactions en- 
hance localization: the higher kinetic energy of localiza- 
tion is compensated by the reduced interaction energy. 
In the extreme limit, a Wigner molecule of electrons is 
formed [isj . It is now tempting to interpret this bal- 
ance between the two localizing effects in the way that 
the interactions affect the anticrossings only indirectly 
via the effective potential. Furthermore, as one cannot 
evidently identify any clear detail of the spectra to result 
from the electron-electron interactions, the possibility of 
FIR spectroscopy to reveal correlation effects is some- 
what questionable. Its ability, however, to probe for the 
effective confinement potential of electrons is clear. This 
is especially true if various polarizations of radiation are 
used. 

To summarize, we have presented FIR spectra of a lat- 
eral two-electron quantum-dot molecule and studied the 
effects of confinement potential and correlations on the 
spectra. We have shown that lowering the symmetry of 
the confinement potential induces changes in the spec- 
tra. These changes are deviations from the two Kohn 



modes, and include the splitting of zero-field resonances, 
anticrossings in w+, and an additional mode above w+ 
in the spin-triplet spectrum. We have studied the effect 
of the electron-electron interactions on these deviations, 
and surprisingly, the interactions drive spectra towards 
those of a high-symmetry QD. All the deviations from 
Kohn modes listed above are reduced as the interactions 
are turned on. This leads to an interesting conclusion: 
the deviations from the Kohn modes result merely from 
the non-parabolic external confinement potential. The 
interactions change the effective potential of the elec- 
trons, but as our results show, the resulting potential 
is more parabolic than the bare external confinement. In 
this way, the resulting FIR spectra are closer to Kohn 
modes than the non-interacting ones. In addition, our 
results show that probing interaction effects by FIR spec- 
troscopy is not straightforward, as clear signatures of in- 
teractions are not present even in the two-electron case 
where the interaction effects should be strongest. On the 
other hand, FIR spectroscopy is able to reveal the non- 
parabolicity of the external confinement potential very 
efficiently, particularly if polarized radiation is used. 
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